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ABSTRACT 



We present high resolution (R=55,000) optical spectra obtained with MIKE 
on the 6.5 m Magellan Clay Telescope as well as Spitzer MIPS photometry and 
IRS low resolution (R~ 60) spectroscopy of the close (14 AU separation) binary, 
HD 101088, a member of the ~12 Myr old southern region of the Lower Centaurus 
Crux (LCC) subgroup of the S cor pius- Centaurus OB association. We find that 
the primary and/or secondary is accreting from a tenuous circumprimary and/or 
circumsecondary disk despite the apparent lack of a massive circumbinary disk. 
We estimate a lower limit to the accretion rate of M > 1 x 10^^ Mq yr~^, which 
our multiple observation epochs show varies over a timescale of months. The 
upper limit on the 70 fim flux allows us to place an upper limit on the mass of 
dust grains smaller than several microns present in a circumbinary disk of 0.16 
Mmoon- We conclude that the classification of disks into either protoplanetary or 
debris disks based on fractional infrared luminosity alone may be misleading. 
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INTRODUCTION 



Most young stars are initially surrounded by optically thick accretion disks (iBeckwith et al. 



19901). Though the frequency of binarity is a function of mass and formati on environm ent, 
one-third of all main sequence stars in the Galactic disk are in binaries (Ladal 120061) . In 
nearby star forming regions, more than half (iGhez et al.lll993l : ISimon et al.lll995l ) of young 
stars have been observed to be members of binary systems. Since the formation of a binary 
is a common outcome of the star formation process, studying the structure and evolution 
of disks in binary systems is an important part of developing a complete understanding of 
the planet formation process. In the case of young binary systems, disks may s urround the 
primary, the secondary, and/or both components (lArtymowicz fc Lubowlll994j ). The sepa- 
ration between stars in a binary system is an important parameter that affects the geometry 
and structure of any disk material surrounding the stars. Optically thick disks around each 
member of young binaries have been observed for stars separated by as little as 14 AU 
( iHartigan fc Kenyonl 120031 ) . The outer radii of disks in close binary systems are truncated 
through gravitational interactions between the binary components, limiting the amount of 
material available for planet formation and accretion onto the stars. 

Optically thick accretion di sks around single stars appear to dissipate within a few 
million years ( iHaisch et al.ll200ll ). Since the disks around each star in close binary systems 
are truncated to smalle r outer radii than disks around single stars, they may disperse faster. 
Bouwman et al.l ( 120061 ) conducted a survey with the Spitzer Space Telescope of the ~8 Myr 
old 7] Chamaeleontis cluster and found that circumstellar disks were detected around 80% 
of single stars yet absent around 80% of the close binary stars. This is suggestive of a 
shorter timescale for disk removal in close binaries although the sample size was small and 
the binaries were not spatially resolved. Additional observations of disks in close binary 
systems are needed to confirm these results. 

Here we present Magellan Inamori Kyocera Echelle (MIKE) R=55,000 optical spec- 
troscopy along with Spitzer Infrared Spectrometer (IRS; Houck et al. 2004) and Multiband 
Imaging Photometer for Spitzer (MIPS; Rieke et al. 2004) observations of the close binary 
star, HD 101088, a member of the ~12 Myr old southern regio n of the Lower Centa urus- 
Crux (LCC) subgroup of the Scorpius- Centaurus OB associati on ( Ide Zeeuw et al.lll999[ ). HD 
101088 consists of a F5 primary star (IHouk fc Cowleyl Il975[ ) and a secondary of unknown 
spectral type according to Hipparcos astrometry. The two components are separated by 
0.15"or 14 AU at a distance of 94 pc (Ivan LeeuwenI 120071 ). Our high resolution optical spec- 
tra reveal broad, spatially unresolved Ha emission from this source which is indicative of 
ongoing stellar accretion. The Spitzer IRS spectrum shows there is very little, if any, hot 
circumstellar dust. The lack of strong emission in the mid- and far- infrared indicates the 
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absence of a cold outer disk. 



OBSERVATIONS AND DATA REDUCTION 



We observed HD 101088 with the MIKE spectrograph (iBernstein et al.ll2003l ) on the 6.5 
m Magellan Clay Telescope at Las Campanas Observatory on 2007 March 11, 2009 April 15, 
and 2009 June 8 (UT). The 0.35" x 5" slit was used in all cases, giving a resolution of 55,000 
at the wavelength of Ha. The spectra were flat-fielded, extracted and wavelen gth calibrated 



using t he MIKE pipe line written by D. Kelson with techniques described in iKelson et al. 



(120001 ) ■ iKelsonI ((20031), and Kelson et al.l fcOOGf ). The exposure times were 24, 60, and 147 
seconds for the 2007 March, 2009 April, and 2009 June observations respectively. The 
resulting spectrum has a signal-to-noise ratio larger than 100 per pixel for wavelengths greater 
than 4500 A. 

Photometry at 24 and 70 /im was obtained using MIPS on Spitzer. The observations 
were made on 2005 April 9 with integration times of one cycle of 3 s at 24 /xm and one cycle of 
10 s at 70 um. The Data A nalysis Tool (DAT), version 2.80, created by the MIPS instrument 
team (jCordon et al.ll2005l ) was used to reduce the data. The MIPS calibration uncertainty 
was taken to be 4% at 24 /im and 7% at 70 /im based on the MIPS handbook. Additional 
MIPS data processing details can be found in C. Chen et al. (2010, in preparation). A 
Spitzer IRS spectrum was obtained on 2007 June 18 with both the Short-Low (5.2-14.0 /tm) 
and the Long-Low (14.0-38.0 /im) modules. The observations were carried out in IRS staring 
mode with no peak-up. The Short-Low observations consisted of 2 cycles of 6 seconds each 
while the Long-Low observations were 2 cycles of 14 seconds. The data were re duced and 
analyzed using the SMART p rogram creat e d by t he IRS team (IHigdon et al.ll2004j ) following 
the procedures described in iFurlan et al.l (120061 ) . with sky subtraction from the opposite 
nod position. We estimate the spectrophotometric uncertainty of the IRS spectrum to be 
approximately 5%. 



RESULTS 



We estimated the age of the system by placing the star on the HR diagram and com- 
paring to pre-main-sequence evolutionary tracks and by considering the age of the ensemble 
of young stars with which HD 1010 88 is associated. Isochr onal age estimates for HD 101088 
range from 2 Myr using the tracks of lPalla fc Stahlerl (120011 ) to 5 Myr based on the lSiess et al. 
(120001 ) tracks. Since observational uncertainties can lead to errors in age estimates for indi- 
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vidual stars, we also co nsidered the average age of th e region of Lower Centaurus Crux where 
HD 101088 is located. iPreibisch &: Mamajekl (120081 ) found an average age for the members 
of the southern portion of Lower Centaurus Crux of 12 Myr so we conservatively estimate 
the actual age of HD 101088 as being between 2 and 12 Myr. 

The projected rotational velocity {v sin i) was determined by fitting rotationally broad- 
ened synthetic spectra to our data. The fitting region was restricted to 4000-7000 A, exclud- 
ing regions with strong telluric contamination, Ha, and the Na D lines. We used Richard 
Gray's spectral synthesis program, SPECTRUM^, along with Kurucz model atmospheres 
of solar metallicity to con apute t he sy nthetic spectra and then broadened them using the 
rotational profile given in iGrayl (Il992[ ). We fitted these synthetic spectra to our data and 
performed a minimization to find the best value for v sin i. The average of the measure- 
ments derived from all the separate observations is listed in Table [U The RVCORRECT and 
FXCOR packages in IRAF were used to measure the radial velocity. Measurements from 
the three different observations were consistent with one another and the average is 17.8 ± 
3.9 km s~^. By combining the position, prope r motion, and parallax of HD 101088 from the 
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revised Hipparcos catalog ( Ivan LeeuwenI 120071 ) with our measured radial velocity, we derive 



a Galactic space velocity (U,V,W) = -5.6 ± 1.9, -21.9 ± 2.8, -6.5 ± 2.1 km s~^, consistent 
with the mean space motions for members of Lower Centaurus Crux (E. Mamajek, private 
communication) . 

The spectra from all observed epochs show broad Ha emission as seen in Figure [T] In 
order to characterize the Ha emission which falls on top of photospheric absorption, we 
subtracted a broadened HD 106444 spectrum, a star with the same spectral type as the 
HD 101088 primary and somewhat smaller v sin i (96 km s~^ vs. 160 km s~^) which we 
also observed using MIKE. The Ha measurements we report are based on the subtracted 
spectra. The Ha emission has an equivalent width of 3.6, 4.3, and 6.2 A in 2007 March, 2009 
April, and 2009 June respectively. We estimate the systematic uncertainty for the equivalent 
width measurements to be 0.4 Awhich is the standard deviation of the measured equivalent 
widths of the Ha absorption in 14 F5 stars we have observed in Sco Cen. The Ha full width 
at 10% of the peak is commonly used to differentiate between accretion and chromospheric 
emission. Ha 10% wi dths > 270 km s~^ are due to accretion independent of spectral type 
dwhite fc Basrill2003h . For HD 101088, the Ha 10% full width is broad but variable over 
our three different observations; 388 ± 2, 380 ± 2, and 429 ± 4 km s"^ in 2007 March, 2009 
April, and 2009 June respectively, suggesting the presence of ongoing accretion. 



We compared the observed Ha profiles with radiative transfer models of magnetospheric 



^http://www.phys.appstate.edu/spectrum/spectruni.htnil 
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accretion ( iMuzerolle et al.ll200ll ). These models are consistent with the h ne emission ob- 



serve d in most classical T Tauri stars, as well as some Herbig Ae stars ( IMuzerolle et al. 



2004J ). and so provide potential constraints on the accretion activity in the HD 101088 sys- 
tem. We calculated models using the mass, radius, and effective temperature of the primary 
star as fixed inputs. The outer radius of the magnetosphere was fixed at the corotation 
radius, which is about 1.5 R* given the v sin i value we measure for the primary. The gas 
temperature, density, and inclination were then varied to find a good fit to the observed 
profile. Figure [T] shows one example fit to the 2007 observation. We are also able to fit the 
profile assuming that the secondary rather than the primary is accreting, using similar pa- 
rameters and accounting for the additional continuum from the unresolved primary. Either 
model has difficulty matching the increased line emission in the 2009 spectrum; the accretion 
geometry may have changed, or an additional source of emission such as a wind may have 
manifested itself. The difficulty in fitting the June 2009 profile is that the models do not 
produce sufficiently strong emission for any reasonable set of parameters. The only way to 
increase the emission would be to increase the size of the magnetosphere, but that is not 
possible since it cannot be larger than the corotation radius. One way around this would be 
to invoke a different geometry from the dipole approximation, such as a pinched or inflated 
field configuration. Note that because of degeneracy in constraining the gas temperature and 
density, we can only put a lower limit to the actual accretion rate, M > 1 x 10"^ Mq yr^^. 
Spatially resolved observations of other diagnostics, such as infrared emission lines or UV 
continuum excess, are needed to obtain better constraints on the accretion rate and its origin 
in the system. 

The Spitzer MIPS photometry gives fluxes of 70.0 ± 2.2 mJy at 24 /im and a 3-0" upper 
limit of 16.5 mJy at 70 /im. The uncertainty for the 24 /zm measurement was computed by 
adding the 0.75 mJy statistical error to the 3% calibration uncertai nty in quadrature. Th e 



primary star in HD 101088 is known to have a spectral type of F5 (IHouk fc Cowley! Il975l ). 
In order to compute the spectral energy distribution for the stars in HD 101088 to look for 
an infrared excess, we need to determine the spectral type of the secondary star. Hipparcos 
gives Hp=6.8±0.07 for the primary and Hp=9.6±0.89 for the secondary. Assuming ari age o f 



5 Myr and knowing that the spectral type of the primary is F5, we used the lSiess et al.l (120001 ) 
evolutionary tracks to translate the magnitude difference between the two stars, including 
uncertainties, to a range for the secondary's spectral type of K0-K5. We computed the 
stellar photospheric flux by normalizing Kurucz model atmospheres for both components of 
the binary to Hipparcos photometric values. We assume HD 101088 has solar metallicity 
and our calculated value for the line-of-sight extinction of A„=0.13. Figure [2] shows the 
combined stellar photospheric flux plotted over fluxes at B and V from Hipparcos and J, H, 
and K from 2MASS. The upper panel shows the SED assuming a K5 secondary while the 
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lower panel shows the SED in the case of a KO secondary overplotted with the IRS spectrum, 
MIPS 24 /im point, and MIPS 70 /xm upper limit. The IRS spectrum appears to have nearly 
the same slope as the stellar photosphere and shows no silicate emission features. In the 
case of a K5 secondary, the MIPS 24 /xm flux is 6.2a larger than the 56.4 mJy flux from 
the combined stellar photospheres. We modeled the infrared excess emission assuming that 
the dust grains are large, consistent with the featureless Spitzer IRS spectrum. However, 
if dust grains smaller than a few microns are present, the lack of spectral features in the 
IRS spectrum requires that the dust grains are composed of material other than silicates, 
perhaps carbon or iron. Because we lack measurements at wavelengths between the 2MASS 
photometric points and the IRS spectrum, acceptable fits can be achieved for a wide range 
of dust temperatures, 500-2000 K. If the secondary is a KO star, then the MIPS 24 /im flux 
is only 1.7a larger than the 66.2 mJy combined stellar photospheric flux and consistent with 
there being no infrared excess. However, with a KO secondary star, the SED overpredicts 
the flux at J, H, and K by 8.4(T, 12.1a, and 4. la respectively. A Kl spectral type is the 
hottest allowable secondary star which agrees with the J, H, and K photometry. Given 
the uncertainty about the nature of the secondary star in the system, we cannot determine 
whether there is an infrared excess due to any warm circumstellar dust. After subtracting 
a F5 template from our HD 101088 spectra, we searched for any leftover absorption lines 
due to the secondary star and found none. This does not necessarily rule out a companion 
since they could be heavily veiled by light from the primary or any accretion excess from the 
secondary if it is accreting. Or, perhaps they are broadened due to rapid rotation. Spatially 
resolved observations are needed to characterize the spectral type of the secondary star. 

The Spitzer MIPS 70 /im upper limit can be used to place a constraint on the presence 
of any cold dust in the system. An upper limit on the mass of circumstellar dust in the 
system contributing to the 70 /im flux can be made assuming the dust is optically thin 
and at a single temperature. In that case, the dust mass is given by M^ust 



( iJura et al.lll995l ). where is the distance to the star and is the dust absorption opacity. 
If we assume a circumbinary disk is present and take the semi-major axis of the binary 
orbit to be half of the Hipparcos separation between binary components, the inner radius of 
the circumbinary disk would be at 15 AU or 2.08 times the semi-major axis of the binary 
orbit, the expe cted location of the inner e dge of a circumbinary disk for a binary with 



circular orbits (lArtymowicz fc Lubowi Il994j ). We estimate a temperature of T = 150K at 
this location by simple radiative balance with the stars, T ~ (L^/lGnr'^aY^^, where a is the 
Stefan-Boltzmann constant. With F^, the upper limit on the 70 /im flu x equal to 16.5 mJy , 



the distance to HD 101088 taken to be 94 pc, and = 3 cm^ g"^ ( jPoUack et al.lll994j ). 
Mdust < 1-2 X lO^^g or 0.16 Mmoon- If dust grains much larger than several microns are 
present in the outer disk, the amount of mass there could be substantially larger. 
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4. DISCUSSION 

The presence of accretion in the absence of a large infrared excess is surprising. Theo- 
retical models have shown that circumstellar disks in close binary systems are expected to 
be tidally truncated to an ou ter radius significantly smaller than disks around single stars 



(lArtymowicz &: Lubowl Il994l ). They give the expected truncation radii for circumprimary 
and circumsecondary disks based on the semi- major axis of the binary orbit. Assuming cir- 
cular orbits, the outer disk radii are 0.46a and 0.2a for circumprimary and circumsecondary 
disks respectively. For HD 101088, this corresponds to outer disk radii of 6.5 AU and 3 AU. 

Since our observations do not spatially resolve the two binary components, we cannot 
decipher which of the two stars has a disk or if both do. However, in either case, the timescale 
to viscously accrete the material fror a such tidally t runca ted disks is very short. We follow 



the viscous accretion calculations in iQuillen et al.l (|2004( ) to compute the time required to 



deplete the disk material. The accretion timescale is given by = a~^{j^YToTb/'^T^^ where a 
is the viscosity parameter, r is the disk outer radius, h is the disk scale height, and Torb is the 
orbital period of material at the outer disk radius. We take a to be 0.0 1, a value typical for 



accretion disks, and h/r is estimated from the parametrization given by lChiang fc Goldreich 
( ]1997n . h/r = 0.17 a /I. For the truncated circumprimary disk, the orbital period at the disk 



truncation radius is 14 years and the time to clear the disk through viscous accretion is 2650 
years. The slightly smaller circumsecondary disk would accrete all its material in only 1800 
years. 

Both timescales are much shorter than the age of the system implying that the circum- 
stellar material must have been replenished, perhaps from a circumbinary disk. There are 
examples of strongly accreting stars which appear to be devoid of material to accrete base d 



on the lack of near- infrared excess emission such as V4046 Sgr (jJensen fc Mathieul 119971 ) , 
however, in contrast to HD 101088, evidence is seen of an outer disk reservoir in the form of 
large excess emission at wavelengths longward of 10 fim. If a reservoir of material exists in 
a circ umbinary disk, material may be transferred to a circumstellar disk around one of the 



stars (lArtymowicz fc Lubowl Il996l ). In this model, circumbinary disks can transfer material 
via pulsed accretion to the circumstellar disks of the individual stars to support accretion 
rates comparable to that found for single stars with full disks. Some observati onal evidence 



for pu lsed accretion which varies as a function of orbital phase has been found. iJensen et al. 



( 120071 ) reported BVRI photometry of the pre-main-sequence spectroscopic binary UZ Tau E 
which varied on a timescale consistent with the period of the binary suggesting that periodic 
accretion can occur from a circumbinary disk. Variability in the Ha emission was not as 
clearly periodic, perhaps due to a lack of data. The variation we observe in the HD 101088 
Ha emission over a timescale of a few months is unlikely to be due to pulsed accretion from 
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a circumbinary disk since the orbital period of the system, ~35 years, is significantly longer. 



It is not clear from our observations that HD 101088 possesses a circumbinary disk from 
which to replenish material around the individual stars. Our Spitzer observations suggest 
there is very little material that could be present in a circumbinary disk. The dust mass esti- 
mate based on the upper limit at 70 /im of < 0.16 Mmnn y. is eight orders of m agnitude smaller 
than the 0.13 Mq circumbinary disk around GG Tau (IDutrey et al.lll994j ). A larger survey 
of circumbinary disk masses has found typical u pper limits of 0.005 Mq for binaries with 
separations less than 100 AU ( jJensen et al.lll996l ). still significantly larger than our derived 
upper limit for HD 101088. Our observations suggest there is no circumbinary disk present 
and we have caught HD 101088 just as it is clearing its remaining circumstellar material. 
Spatially resolved observations to determine which star is accreting and submillimeter data 
to probe any outer disk material present are needed to further elucidate the properties of 
this interesting system. 



5. CONCLUSIONS 

We have obtained high resolution optical spectroscopy with MIKE on the Magellan Clay 
telescope as well as Spitzer IRS spectroscopy and MIPS 24,70 fj,m photometry of the close 
binary, HD 101088. 

1. Broad Ha emission is present in our spectra and reveals ongoing stellar accretion 
and the presence of circumstellar gas. We derive a lower limit on the accretion rate of M 

> 1 X IQ-^ Mq yr-\ 

2. The truncated disks in such a close binary have viscous accretion lifetimes much 
shorter than the age of the system, suggesting that some source of replenishment must be 
present to maintain the ongoing accretion. 

3. The upper limit at 70 /im leads to a constraint on the dust mass in a circumbinary 
disk if it is present of < 0.16 Mmoon- 

4. HD 101088 would be classified as having a debris disk based on its small fractional 
infrared luminosity, 7.0x10"^. The presence of ongoing accretion shows that the classification 
of disks into either protoplanetary or debris disks based on fractional infrared luminosity 
alone may be misleading. 

5. We find that the primary and/or secondary is accreting from a tenuous circumprimary 
and/or circumsecondary disk despite the apparent lack of a massive circumbinary disk. This 
unique situation merits further study. 
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Fig. 1. — Three Ha line profiles obtained during different epochs of observation overplotted 
with an example accretion model line profile (dashed line) for the March 2007 observation. 
The parameters for the model line profile are M = 1 x 10^^ Mq yr~^, Tma^=10,000 K, i=10°, 
and Rmag=l-3-1.5 R*. Due to degeneracy in constraining the gas temperature and density, 
we can only put a lower limit to the actual accretion rate, M > 1 x 10~^ Mq yr~^. The raw 
March 2007 spectrum around Ha before subtracting the template spectrum to remove the 
underlying photospheric absorption is shown offset from the three line profiles. 

Fig. 2. — Possible SEDs for the range of secondary spectral types allowed by Hipparcos 
observations and uncertainties. The top panel assumes a K5 secondary and shows the total 
contribution from both components of the HD 101088 binary plotted over B and V fluxes 
from the Tycho-2 Catalogue and 2MASS J,H, and K fluxes. The error bars on the B,V,J,H, 
and K values are smaller than the plotted symbols. The Hipparcos Hp magnitude of the 
secondary star is plotted along with its error bar. The dashed lines show the individual 
contributions from the two stars, a F5 primary and K5 secondary. The Spitzer IRS spectrum 
is shown along with the MIPS 24 /xm photometric point and MIPS 70 fim 3a upper limit. 
The inset shows a region around the MIPS 24 /xm point with the la error bar overplotted. 
The combined photospheric flux in this case is 6.2a below the MIPS 24 fim point. The 
bottom panel shows the same with a KO secondary. In this case, the combined photospheric 
flux is only 1.7cr below the MIPS 24 /xm point but above the 2MASS J,H,K points by 8.4cr, 
12.1(7, and 4.1(7 respectively. 



This preprint was prepared with the A AS macros v5.2. 
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Table 1. Properties of HD 101088 



Property 


Value 


Reference 


Primary Spectral Type 


F5IV 


1 


Secondary Spectral Type 


KO-KS'^ 


2 


Distance (pc) 


9411 


3 


Age (Myr) 


2-12*^ 




RV (km s-i) 


17.8 ± 3.9 


2 


V sin i (km s~^) 


160 ± 4 


2 


Ha EW (A) 


3.6/4.3/6.2^-^ 


2 


Ha 10% FW (km s"^) 


388/380/429^ ^ 


2 


L/r/L* 


7.0x10"^ 


2 



'^Spectral type range based on Hipparcos magnitude 
and uncertainties. 

'^2-5 Myr from HR diagram fits, 12 Myr average age of 
southern LCC 

'^Ilmar07/15apr09/08jun09 

■^Ha EW Systematic Uncertainty: 0.4 A 

^Ha 10% FW Uncertainties: 2/2/4 km s^^ 



Referen ces. — (1) iHouk fc CowlevI (119751 ). (2) this 
work. (3) van LeeuwenI ( 2007 ). 



